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ABSTRACT 


The  present  concept  of  automatic  throttle  control,  as  employed  in 
Navy  carrier-based  aircraft,  was  investigated.   The  aircraft  chosen  for 
study  was  the  A-7E.   The  powerplant  was  the  TF41-A-2,  a  turbofan  engine 
with  a  relatively  slow  throttle  response  in  the  approach  power  range. 

The  effects  of  additional  inputs  to  the  approach  power  compensator 
were  evaluated.   It  was  shown  that  a  considerable  increase  in  performance 

could  be  achieved  through  the  incorporation  of  longitudinal  feedback. 

In  addition,  the  limitations  imposed  on  performance  by  large  engine  lags 

were  found  to  be  much  less  severe  for  systems  with  longitudinal  feedback. 

The  modifications  suggested  require  a  redesign  of  the  approach  power 

compensator  system  currently  in  use  by  the  Navy. 
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I.   INTRODUCTION 

A.   APCS  PERFORMANCE  OBJECTIVES 

The  primary  purpose  of  an  approach  power  compensator  system  (APCS) 
in  Navy  carrier-based  aircraft  is  to  provide  the  thrust  necessary  to 
maintain  the  correct  airspeed  during  carrier  approach  and  landing.   The 
APCS  can  be  used  either  in  an  automatic  carrier  landing  system  (ACLS) 
or  in  a  piloted,  manual  approach.   Although  airspeed  systems  exist 
[Ref.  l] ,  the  systems  utilized  in  Navy  aircraft  attempt  to  maintain  a 
constant  angle  of  attack  (AOA) .   Supplementary  inputs,  such  as  normal 
acceleration  and  elevator  crossfeed,  are  usually  employed  in  an  effort 
to  maintain  the  "on-speed"  AOA.   Since  the  on-speed  value  of  AOA  does 
not  vary  from  approach  to  approach,  there  is  no  indexing  requirement. 
In  an  airspeed  system,  however,  the  reference  airspeed  is  a  function 
of  the  aircraft  configuration  in  the  landing  approach.   As  configuration 
variables,  such  as  weight,  vary  from  approach  to  approach,  so  does  the 
reference  airspeed.   This  is  an  important  factor  in  favor  of  an  AOA 
system.   The  pilot  workload  which  airspeed  indexing  would  add  to  an 
already  demanding  task  probably  would  not  be  tolerated,  and  the  designer 
of  an  approach  power  compensator  system  should  bear  this  in  mind. 

Military  specifications  reflect  the  objective  of  maintaining  a 
constant  AOA.   Reference  2  states  that  the  output  of  the  APCS  shall  be 
proportional  to  error  in  AOA,  change  in  normal  acceleration,  integral 
of  AOA  error,  and  elevator  position.   As  is  noted  in  Ref.  1,  it  appears 
that  the  APCS  design  concept  is  overly  restricted  by  military  specifi- 
cations.  However,  these  specifications  do  point  out  the  current  Navy 
APCS  concept. 
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B.  AIRCRAFT 

The  aircraft  chosen  for  study  was  the  A-7E,  a  light-attack,  carrier- 
based  aircraft  of  approximately  20,000  pounds  empty  gross  weight.   The 
A-7E  was  chosen  because  it  was  one  of  two  aircraft  included  in  a  previous 
study  of  Navy  approach  power  compensator  problems  and  requirements 
[Ref.  1]  and  because  the  author  has  had  personal  experience  as  an  A-7A 
pilot. 

The  A-7E  powerplant  is  a  turbofan  engine  whose  throttle  response 
is  characterized  by  a  relatively  long  and  nonlinear  time  constant  in 
the  approach  thrust  range.   The  APCS  is  an  AOA  system  with  inputs  of 
AOA,  integral  of  AOA,  normal  acceleration,  and  unit  horizontal  tail 
(UHT)  movement.   A  detailed  description  of  the  system  is  contained  in 
Ref.  3. 

C.  PREVIOUS  STUDIES 

Systems  Technology,  Incorporated  conducted  a  study  of  Navy  approach 
power  compensator  problems  and  requirements  [Ref.  l].  The  following  is 
a  summary  of  pertinent  conclusions  and  recommendations: 

1.  A  fundamental  conceptual  difficulty  with  current  APCS's  is 
the  use  of  AOA  feedback  to  constrain  both  AOA  and  airspeed,  as 
the  two  are  only  in  phase  at  low  frequencies.   Thus,  it  is 
possible  to  control  only  one  of  these  or  a  linear  combination 
of  the  two  with  a  single  control  input.   As  a  consequence,  APCS 
design  involves  a  compromise  between  these  constraints. 

2.  Airframe  characteristics  restrict  the  ability  to  achieve 
satisfactory  performance.   Engine  response  time  is  a  fundamental 
problem  area,  as  the  APCS  cannot  augment  the  dynamics  without 
excessively  overdriving  the  engine.   It  was  recommended  that 
the  effect  of  engine  lag  on  APCS  performance  be  investigated 

in  detail  with  a  view  toward  specifying  necessary  characteristics. 


3.  Gust  response  within  the  current  APCS  concept  will,  in 
general,  be  poor.   Gust  proofing  against  longitudinal  gusts 
with  an  AOA  system  is  difficult.   Gust  proofing  against 
vertical  gusts  has  minimal  effect,  as  the  dominant  response 
is  determined  by  the  short  period  characteristics. 

4.  It  was  recommended  that  the  effects  of  additional 
feedbacks  to  the  APCS  be  analyzed  in  relation  to  gust 
response  and  to  optimum  performance  for  separate  con- 
figurations of  fully  automatic  and  manually  controlled 
systems. 

The  Naval  Air  Test  Center  conducted  an  evaluation  of  the  APCS  in 

the  A-7E  [Ref.  3],   The  emphasis  was  on  optimizing  the  performance  of 

the  APCS  for  both  manual  and  ACLS  approaches  with  minimal  material  and 

design  modifications.   The  recommended  modifications  were  basically 

changes  in  feedback  gains  and  the  incorporation  of  a  dual  time  constant 

in  the  UHT  crossfeed  circuit.   A  shorter  time  constant  was  used  for 

nose  up  corrections  than  for  nose  down  corrections,  with  the  net  result 

that  the  UHT  input  was  more  effective  for  nose  up  than  for  nose  down 

attitude  changes. 

D.   OBJECTIVES 

The  purpose  of  this  research  was  to  improve  the  performance  of  the 
APCS  of  the  A-7E  aircraft  as  used  in  the  ACLS.   Specifically,  the 
effects  of  additional  inputs  to  the  APCS  computer  were  evaluated  with  a 
goal  of  improving  the  current  concept  of  APCS  design.   In  addition,  the 
effect  of  engine  response  time  on  APCS  performance  was  determined. 
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II.  METHOD  OF  ANALYSIS 

A.   DESCRIPTION  OF  THE  MODEL 

The  fully  automatic  carrier  landing  mode  was  selected  for  study  in 
this  analysis  of  the  APCS.   The  manual  mode  with  the  pilot  in  the  loop 
was  not  considered.   Thus,  the  total  closed  loop  system  consisted  of 
the  airframe,  the  engine,  the  approach  power  compensator,  the  longitu- 
dinal automatic  flight  control  system  (AFCS) ,  and  the  automatic  carrier 
landing  system. 

The  aircraft  equations  of  motion  were  linearized  about  the  steady 
state  approach  conditions  in  accordance  with  standard  small  perturbation 
theory,  as  described  by  Etkin  [Ref.  4],   Lateral-directional  dynamics 
were  not  included.   The  resultant  equations  are  shown  in  Figure  1.   Air- 
craft stability  derivatives  and  other  parameter  values  were  taken  from 
Ref.  1  and  are  listed  in  Table  I.   Force  and  moment  stability  deriva- 
tives are  normalized  with  respect  to  mass  and  moment  of  inertia,  respec- 
tively.  Hence,  mass  and  moment  of  inertia  do  not  appear  explicitly  in 
the  equations  of  motion,  as  they  do  in  Etkin' s  notation. 

The  A-7E  engine  is  the  TF41-A-2  turbofan.   Reference  5  contains  an  ■ 
analysis  of  the  thrust-power  lever  relationship  for  a  range  of  operating 
conditions.   The  nominal  approach  thrust  for  the  operating  conditions 
as  shown  in  Table  I  is  3000  pounds.   In  this  range  of  thrust,  the  engine 
response  time  constant  is  a  nonlinear  function  of  thrust.   The  relation- 
ship between  the  engine  time  constant,  T  ,  and  thrust,  and  the  relation- 
ship between  thrust  and  power  lever  angle  (6 _.)  for  the  given  operating 
conditions  are  taken  from  Ref.  1  and  are  shown  in  Figure  2.   An  average 
value  of  275  pounds  thrust  per  degree  6pTA  was  used  in  the  model.   Thus, 
the  variable  engine  lag  was  the  only  nonlinearity  in  the  engine  model. 
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Figure  3  is  a  schematic  of  the  APCS  loop  closures.   In  small  per- 
turbation theory,  AOA  is  defined  as  w/U~,  in  radians.   Thus,  the  inputs 
to  the  APCS  computer  are  proportional  plus  integral  w,  a1  (vertical 
acceleration  corrected  for  accelerometer  location),  and  a  filtered  UHT 
feed f orward  term. 

The  APCS  gains  are  a  representative  set  listed  as  "PAX"  gains  in 

Kef.  1  and  do  not  correspond  exactly  to  those  in  use  in  the  current  A-7E 

APCS  configuration.   In  addition,  several  simplifications  were  made  to 

the  APCS  model.   A  complex  pole  at  ten  rad/sec  in  the  throttle  actuator 

response  was  neglected.   This  frequency  is  well  beyond  the  range  in 

which  the  APCS  is  effective.   Throttle  linkage  hysteresis  of  0.5  degrees 

was  ignored.   The  implementation  of  a  dual  UHT  input  time  constant,  as 

described  in  the  introduction,  was  not  incorporated  into  the  model.   A 

gain  adjust  bias  exists  in  the  APCS  which  changes  all  APCS  gains  by  the 

same  amount  to  compensate  for  ambient  air  temperature  effects.   The 

standard  day  value  of  1,0  was  used.,   K   ,  the  vertical  acceleration  gain, 

n 
z 

has  two  values,  the  smaller  value  being  for  load  factors  in  excess  of 
1.1  g's.   Only  the  larger  value  was  incorporated  into  the  model. 

Both  the  longitudinal  AFCS  and  the  ACLS  models  were  taken  from 
Ref.  6.   As  above,  the  AFCS  and  ACLS  gains  were  representative  but  were 
not  necessarily  identical  to  those  currently  in  use.   A  complex  pole 
and  a  first  order  lag  at  20  rad/sec  in  the  AFCS  response  was  neglected. 
The  AFCS  output  is  a  UHT  deflection,  which  is  a  function  of  attitude 
error  (9-G  ),  pitch  rate,  and  normal  acceleration.   The  automatic  carrier 
landing  system  was  represented  by  the  A-7E  SPN-42  longitudinal  control 
equation.   The  output,  a  ship-to-aircraft  pitch  command,  is  a  function 
of  the  aircraft  altitude  error,  Z  ,  measured  perpendicular  to  the  ideal 
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glide  slope.   In  practice,  the  aircraft  slant  range  and  elevation  are 
measured  by  the  SPN-42  radar,  corrected  for  ship  motion  and  radar 
location,  transformed  into  cartesian  coordinates,  and  utilized  to 
command  an  aircraft  pitch  angle. 

The  AFCS  and  ACLS  longitudinal  control  equations  are  shown  in 
Figure  4.   APCS,  AFCS,  and  ACLS  gains  are  listed  in  Table  II. 

B.   DIGITAL  SIMULATION 

1.  Continuous  System  Modeling  Program 

The  closed  loop  system  was  simulated  on  the  IBM  360  digital  computer 
utilizing  a  digital  simulation  language,  CSMP.   CSMP  is  an  acronym  for 
Continuous  System  Modeling  Program.   The  program  is  augmented  by  basic 
FORTRAN  and  provides  a  set  of  functional  blocks  which  simulate  such 
analog  components  as  integrators,  relays,  and  function  generators.   A 
detailed  description  of  the  program  is  contained  in  Ref,  7.   A  descrip- 
tion of  CSMP  functions  used  in  the.  simulation  is  provided  in  Table  III. 

2.  CSMP  Program  Components 

A  sample  program  and  output  are  contained  in  Appendix  A.   The 
program  consists  of  five  major  sections:   the  aircraft  equations  of 
motion  and  engine  model,  the  APCS  model,  the  AFCS  model,  the  ACLS  model, 
and  a  gust  input  model. 

a.   Aircraft  Equations  of  Motion  and  Engine  Model 

The  equations  of  motion  as  shown  in  Figure  1  were 
put  into  state  format: 

{X}  =   [A]  {X}   +   [B]  {R} 

(1) 
{Y}  =   [C]  {X}   +   [D]  {R}   +   [E]  {Y} 
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where: 


{X}   - 


W 

e 
q 


{y}  = 


Cr) 


_/AT 


a' 
z 


w 


g 


(2) 


The  A,  B,  C,  D,  and  E  matrices  are  listed  in  Table  IV. 

The  engine  modeling  consisted  of  two  parts.   As  mentioned  in 
section  II. A.  ,  a  relationship  of  275  pounds  thrust  per  degree  6.^..  was 
assumed.   The  engine  time  constant,  as  shown  in  Figure  2,  was  approx- 
imated by  a  linear  function  generator  of  fifteen  unequally  spaced  points. 
The  thrust  operating  range  modeled  was  ±2000  pounds  about  the  nominal 
operating  value  of  3000  pounds. 

b.   APCS  Model 

The  APCS  equations  were  obtained  from  Figure  3.   The 

throttle  actuator,  power  lever,  and  engine  gains  were  combined  into  a 

single  parameter.   As  a  result,  the  variable  labeled  "PLA"  in  the  CSMP 

program  is  not  the  actual  power  lever  angle.   Gains  and  other  constants 

in  the  feedback  loops  were  also  consolidated  wherever  possible.   AOA 

and  integral  of  AOA  were  kept  distinct.   For  ease  of  programming,  the 

6  /6   transfer  function  was  written  as 
e   e 


(s)   = 


K6  T  s 
e  w 

T  s+1 
w 


s  +  1/T 


w 


V 


wi 

6 


Wl 


(3) 


The  UHT  input  equations  in  the  simulation  are  in  the  form  of  the  right 
side  of  Equation  (3). 
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c.  AFCS  Model 

The  AFCS  equations  were  obtained  from  the  equation  given 
in  Figure  4  with  numerical  values  substituted  for  the  various  gains. 

d.  ACLS  Model 

The  altitude  error,  Z  ,  was  defined  as  the  negative  of  h, 
as  the  two  are  measured  in  opposite  directions  perpendicular  to  the 
ideal  glideslope.   The  6  /Z   transfer  function,  as  shown  in  Figure  4, 
was  rearranged  for  programming  as  follows: 


e 

F(s) 

e 


K  K 
c  x 

K  s+1 
P 


3.9 


+  1 


M 


+ 


3.5 


+  1 


t.K 
-^*  + 


+ 


R  s 
x 


w 


+  ts  +  \ 


(4) 


=  K 


where 


s  +a  s+an 


s(s+l/K  ) (s  +7.14s+25.0) 


K  = 


K  K  (97. 5R  t  +25.0) 
c  x  x  r J 

3.9K 


3.9+t.K 

L  X 

ll     3.9R  t  +1.0 
x  r 


3.9t.K 
_  1  x 

l0     3.9R  t  +1.0 
x  r 


Finally 


■<s) 


K  .    S  +aiS+a0 

s  +b_s  +b  s+b 


Z    Z. 
e    i 


(5) 


where 


K 


+  7.14      b.  =  —^     +  25.0 
1      K 


25.0 
K 
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The  SPN-42  equations  in  the  simulation  are  in  the  form  of  the  right 
side  of  Equation  (5). 

e.  Gust  Model 

To  simulate  a  random  gust  field,  the  output  of  a  random 
number  generator  with  a  uniform  amplitude  distribution  was  passed  . 
through  a  second-order  filter  with  a  break  frequency  of  1.0  rad/sec  and 
a  damping  ratio  of  0.707.   The  theoretical  derivation  of  the  spectral 
properties  of  the  resultant  signal  is  described  in  Appendix  B.   Vertical 
and  horizontal  RMS  gust  velocities  of  5.0  ft/sec  were  used  in  the  simu- 
lation.  To  insure  identical  gust  signals  for  each  run,  a  constant  step 
size  was  used  in  the  numerical  integrations  for  all  runs. 

f.  RMS  Velocity  Error  Criteria 

An  RMS  velocity  error  for  u  and  w  was  used  as  a  measure 
of  APCS  effectiveness.   It  was  defined  as: 


^Mi 


T 
i  T[u(t)2  +w(t)2]dt 


(6) 


where  T  is  the  run  length. 

g.   Digital  Simulation  Error  Analysis 

A  fourth-order  Runge-Kutta  numerical  integration  scheme 
with  a  fixed  step  size  was  used  in  the  simulation.   Figure  5  is  a  plot 
of  normal  acceleration  after  0.5  seconds  versus  step  size  for  an  initial 
altitude  error.   Plots  of  other  system  variableSj  after  0.5  seconds  and 
at  other  values  of  time,  show  the  same  trend.   A  step  size  of  0.05 
seconds  was  chosen  for  the  simulation. 
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C.   THEORETICAL  ANALYSIS 

The  theoretical  analysis  was  primarily  a  root  locus  study  of 
various  APCS  configurations.   To  reduce  the  multiple  input,  multiple 
output  control  system  to  a  form  adaptable  to  the  analysis,  a  multiple 
loop  analysis  technique,  documented  in  Refc  8,  was  employed.   An  expla- 
nation of  the  method,  including  an  example  which  covers  all  applications 
used  in  this  study,  is  contained  in  Appendix  C.   Table  V  is  a  listing 
of  the  open  loop  aircraft  numerators,  cross  product  numerators,  and  the 
open  loop  denominator  used  in  the  analysis. 
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III.  PROCEDURE 

A.   DIGITAL  SIMULATION  RUNS 

.A  standard  set  of  runs  for  the  digital  simulation  was  utilized  for 

each  configuration  of  the  APCS.   The  APCS  configuration  shown  in  Figure 

3  is  referred  to  as  the  standard  configuration.   The  length  of  the 

simulation  runs  was  chosen  to  approximate  the  final  approach  phase  of 

the  carrier  landing.   Under  SPN-42  control,  two  of  the  ACLS  equation 

gains,  R  and  K  ,  are  functions  of  range  for  a  range  greater  than 

6000  ft.   For  ranges  less  than  6000  ft.  the  values  of  R  and  K  are 
°  xx 

constant  at  the  values  shown  in  Table  II»   The  ACLS  equation  employed 
was  valid  for  range  less  than  6000  ft.  to  the  point  at  which  deck 
motion  compensation  is  introduced.   A  run  time  of  30  seconds  was  chosen 
to  approximate  this  phase  of  the  approach.   A  set  of  six  standard 
simulation  runs  was  chosen  to  simulate  system  response  to  initial  high 
and  low  airspeed  and  altitude  errors  and  to  vertical  and  longitudinal 
gust  inputs.   The  set  of  conditions  for  the  six  runs  were: 
Run    (1)    u(0)   =    5  ft/sec,  w(0)  -  -5  ft/sec 

(2)  u(0)   =  -  5  ft/sec,  w(0)  =  5  ft/sec 

(3)  Ze(0)  =   10  ft 

(4)  Ze(0)  -  -10  ft 

(5)  u     =    5  ft/sec 

o 

(6)  w     =    5  ft/sec 

Initial  conditions  and  RMS  gust  velocities  not  specifically  indicated 
in  a  run  are  zero. 
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B.   ANALYSIS  OF  THE  STANDARD  APCS  CONFIGURATION 

To  apply  the  analysis  technique  described  in  Section  II.  C,  the  APCS 

loop  closures  as  shown  in  Figure  3  were  reduced  to  an  equivalent  set 

that  conformed  to  the  block  diagram  format  shown  in  Figure  C.l.   The 

equivalent  transfer  functions  were  a  consolidation  of  all  the  elements 

in  each  loop  closure,  including  the  thrust  to  power  lever  terms.   The 

resultant  diagram  is  shown  in  Figure  6.   G  represents  the  proportional 

w 

plus  integral  AOA  feedback  transfer  function  and  is  given  by 


350.7 


w  u_     IT  s+1    s  /  T  (s+l/T  ) 

e      e 

)|#  +  K   |(s  + 


T     aT/V    K  +K  T 

o>  1/  \    <y     a     <y 


\Jn   T  s(s+l/T  )  (s+l/T  ) 
0  e      cr      e 

Numerically 


(7) 


C 


1021. 5(s+0. 197) 

w        T  s(s+2.857)(s+l/T  ) 
e  e 


G  represents  the  normal  acceleration  feedback  transfer  function.   The 
z 

effect  of  the  6.7  ft.  offset  of  the  accelerometer  on  APCS  performance 
has  been  shown  to  be  minor  in  the  low  frequency  range  [Ref.  1],   Thus, 
the  offset  was  neglected  in  the  equivalent  system,  simplifying  the 

analysis. 

K   -350.7 
n 

G 


z        32.2-T  (s+l/T  ) (s+l/T  ) 
e      a       e 


(8) 

327.2 


T  (s+1)  (s+l/T  ) 
e  e 
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G.   represents  the  UHT  crossfeed  transfer  function  and  is  given  by 
6 

e 

Kc  -350.7s 
o 

G 


6        T  (s+l/T  ) (s+l/T  ) 
e        e      vr    v     e 


(9) 


5(L?_235s 


T  (s+0. 435)  (s+l/T  ) 
e  e 


The  theoretical  analysis  was  confined  to  the  APCS  loop  closures  only; 

AFCS  and  ACLS  loop  closures  were  not  considered.   Justification  for  this 

procedure  is  given  in  Section  IV. B. 

The  closed  loop  transfer  functions  u/6  ,  w/6  ,  and  u/u   for  the 

e'     e'        g 

standard  APCS  configuration  were  formed.   The  engine  time  constant,  T  , 
was  kept  as  a  parameter,  and  root  locus  plots  of  the  zeroes  of  the 
closed  loop  denominator  and  numerators  were  plotted  to  show  the  effects 
of  engine  lag. 

Since  there  are  two  feedback  paths,  as  shown  in  Figure  6,  the 
closed  loop  denominator  is  defined  as  A1  ' .   It  was  derived  by  applying 
equation  C .5. 

a"  -  a  +  VIt.  +  gznIt  <10> 

a 
N       was   expanded   as 

a  (sw-U0s9) 

V  =   NAT  "   SNAT   "  VV  (U) 


Equation  11  was  substituted  into  Equation  10  to  yield 

A"  =a+g„nIt  +  sgz(nIt-VaV  (12) 
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The  expression  was  cleared  over  a  common  denominator,  and  the  numerator 
expressed  in  two  terms,  one  a  function  of  T  ,  to  give 


U  +  1/T  V 
e 


s(s+l)(s+2.857)(s+l/T  ) 


(13) 


where 


U  =   s2(s+l)(s+2.857)  A 


V  ='  s7+4.913s6+9.241s5+10.91s4+6.153s3+1.240s2 


+0.1810s-0. 007699 


The  zeroes  of  A"  were  plotted  as  a  function  of  1/T  by  forming  an 

artificial  transfer  function,  G,  given  by 

1/T  V 
G  -  ~f-  (14) 


G  represents  the  open  loop  transfer  function  in  a  unity  feedback,  single 
loop  system.   This  form  was  adaptable  to  a  documented  root  locus  digital 
computer  program  [Ref.  9].   The  closed  loop  characteristic  equation  for 
the  artificial  system  was 

A'   =  U  +  1/T  V  (15) 

1/T  equal  to  zero  represents  an  infinite  engine  lag,  hence,  a  no  thrust 
condition.   Thus,  the  root  loci  originate  at  the  open  loop  phugoid, 
short  period,  and  control  equation  roots.   The  zeroes  of  the  artificial 
closure  represent  the  roots  that  would  be  obtained  for  the  given  APCS 
feedback  gains  if  there  were  no  engine  delay.   Figure  7  is  a  plot  of 
the  root  loci.   Only  the  loci  above  and  on  the  real  axis  are  shown.   An 
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engine  delay  of  1.8  seconds  occurs  at  the  nominal  operating  thrust  of 

3000  pounds.   This  value  and  the  values  at  thrust  levels  of  ±2000  pounds 

about  the  nominal  value  are  indicated  on  Figure  7. 

Closed  loop  numerator  zeroes  were  plotted  for  N, 

c 

std  std  std 
where  the  subscript  "std"  refers  to  the  loop  closures  for  the  standard 
configuration. 


u~l    w"l     ,  u~\ 


std 


a 
u       u   z      u  w        u 
=  Nc   +  G  Ne   *  +  G  N.  im  +  G.  Nim 

6     z  6  AT    w  6  AT    6  AT 
e       e         e       e 


U  +  1/T  V 

e __ 

s(s+l)  (s+2.857)  (s+0.435)  (s+l/T  ) 


where 


U  -   s2(s+l) (s+2.857) (s+0. 435)  n" 

o 
e 

V  =   -65.43s7-318.6s6-542.8s5-1024.6s4-718.8s3 

-87. 79s2+0. 2346+1. 707 


Figure  8  is  a  plot  of  the  zeroes  of  Ne  | for  1/T  variation. 
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U  =  s(s+1)(s+0.435)nY 

o 
e 

V  -  -2.155s5-511.6s4-782.7s3-343.6s2-68.90s-8.753 


Figure  9  is  a  plot  of  the  zeroes  of  Nc  J  for  1/T  variation, 


of  N™~|  for  1, 
std 


a 
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g       g         g 


U  +  1/T  V 
e 
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g 
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Figure  10  is  a  plot  of  the  zeroes  of  N_  |  for  1/T  variation. 


■oes  of  N   for  1, 


std 

C.   INVESTIGATION  OF  ADDITIONAL  APC5  INPUTS 

In  studying  the  effects  of  additional  APCS  feedback  variables,  T 
was  fixed  at  the  nominal  value  of  1.8  seconds.   The  closed  loop  denomina- 
tor and  numerators  evaluated  at  this  value  are  shown  in  Table  VI.   The 
additional  feedback  variables  which  were  considered  were  combinations  of 
proportional,  integral  and  derivative  of  longitudinal  velocity,  u.   As 

a  result,  the  closed  loop  u/u   and  u/6   transfer  function  zeroes  were 

g        e 

not  affected.   Subsequent  to  each  root  locus  analysis,  a  set  of  digital 
simulation  runs,  as  described  in  Section  III. A,  were  completed  with  the 
altered  configuration.   A  summary  of  the  RMS  velocity  errors  is  contained 
in  Table  VII. 
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Longitudinal  acceleration  feedback  was  first  considered.   The 
equivalent  transfer  function  for  u  feedback,  defined  as  G*  ,  is 


350. 7K- 

c.   =  i_ 

u     T  (s+l/T  ) 


With  T  equal  to  1.8  seconds,  this  becomes 
194. 8K- 


The  closed  loop  denominator,  defined  as  A''1,  was  formed  by  modifying 
A11  as  follows: 

A'"   =  A'  '   +  G-Nf"  (20) 

u  AT 

Figure    11    is   a  root    locus   plot   of   the   zeroes   of  A1  ' '    as  K*   was  varied. 

u 

In  a  like  manner,  longitudinal  velocity  feedback  was  considered. 

194. 8K 
u     s+0.56 

The  closed  loop  denominator  was  again  defined  as  A'  ' ' . 


A"'   =  A"   +  G  N?  •  (22) 

u  AT  v  ' 

Figure   12   is   a   plot  of   the   zeroes   of  A"  '    as  K     was  varied. 

u 

A  combination  of  proportional  plus  integral  longitudinal  velocity 
was  considered  next. 

G"   -  l9A'l^        '    K  <1+K  /s)  (23) 

u     s+0.56      up 


Values  of  K   from  0.1  to  0.5  were  considered.   Figure  13  is  a  plot  of  the 
P 

zeroes  of  A'"  as  K  was  varied  with  K   equal  to  0.5. 

u  P 
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The  final  configuration  considered  was  proportional  plus  derivative 
of  longitudinal  velocity. 


G_  =  l9/*'lr      *   K-(s+K  )  (24) 

u     s+0.56      u    q 


Figure  14  is  a  plot  of  the  zeroes  of  A'  ' '  as  K*  was  varied  with  values 

of  K  equal  to  0.2,  0.5,  and  1.0.   The  phugoid  roots  are  indicated  for 

all  three  values  of  K  .   The  closed  loop  numerator  for  the  w/6   transfer 

q  e 

function  was  formed. 


<   =  rff  +  G_nY  U  (25) 

6e    u  6e£T 


rev 

The  subscript  "rev"  refers  to  revised  configuration.   The  resultant 
root  loci  for  K  equal  to  0.5  are  shown  in  Figure  15. 

D.   REVISED  APCS  CONFIGURATION 

Based  on  the  results  of  the  investigation  of  various  inputs,  a 
final  APCS  configuration  was  chosen.   As  compared  to  the  standard 
system,  this  configuration  had  additional  inputs  of  proportional  plus 
derivative  of  longitudinal  velocity  with  specific  gains  of  0.5  and 

2.0  for  K  and  K* ,  respectively.   The  revised  equivalent  APCS  block 

q      u '     r       ■>  "» 

diagram  is  shown  in  Figure  16. 

With  the  gains  fixed  at  the  values  given  above,  root  locus  plots 


1  „  t, 


for  the  zeroes  of  A"'  and  N    for  1/T  variation  were  plotted  and 

6e 


rev 
are  shown  in  Figures  17  and  18,  respectively. 

Straight-line  approximations  to  Bode  plots  for  u/5   w/6  ,  and 

e'     e 

u/u  were  constructed  for  comparing  the  standard  and  revised  APCS 

O 
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configurations.   The  transfer  functions  are  listed  in  Table  VIII,  and 
the  magnitude  portion  of  the  Bode  plots  are  shown  in  Figures  19  through 
21. 
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IV.  DISCUSSION  AND  RESULTS 

A.   EFFECTS  OF  ENGINE  RESPONSE 

Strictly  speaking,  the  incremental  root  locus  technique  showing  how 
the  poles  and  zeroes  of  the  APCS-equipped  aircraft  vary  with  the  non- 
linear engine  time  constant  is  not  valid.   The  use  of  root  locus  analyses 
is  valid  only  for  linear  systems.   However,  an  incremental  root  locus 
technique  ,has  been  used  for  stability  analyses  of  systems  with  single- 
valued  continuous  nonlinearities  [Ref.  10],   The  system  response  is 
approximately  defined  by  the  instantaneous  location  of  the  roots,  here 
a  function  of  the  instantaneous  thrust. 

The  initial  conditions  and  the  gust  inputs  used  in  the  simulation 
were  considered  to  be  representative  of  those  that  an  APCS/aircraft 
system  would  be  expected  to  encounter.  Visual  inspection  of  system 
response  curves  did  not  reveal  -the  thrust  response  nonlinearity ;  i.e., 
the  curves  were  approximately  sinusoidal  in  shape,  indicating  that  the 
system  performance  could  be  evaluated  by  considering  a  fixed  value  of 
T  equal  to  1.8  seconds.   The  minor  effect  of  the  nonlinearity  was  due 
to  the  low  variance  of  AT  in  relation  to  the  size  of  the  thrust 
operating  range  modeled.   The  dependence  of  the  response  on  T  was 
greater  for  the  revised  system, as  can  be  seen  by  comparing  Figures  7  and 
9  with  Figures  17  and  18.   This  was  a  direct  result  of  the  increased 
capability  to  alter  APCS  performance  with  the  longitudinal  feedback. 

The  analysis  did  not  provide  insight  concerning  the  specific 
effects  of  a  large  engine  lag  on  APCS  performance.   However,  the 
limitations  imposed  on  performance  by  a  large  engine  lag  were  shown  to 
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be  less  severe  for  a  system  incorporating  longitudinal  feedback.   As 
reported  in  Ref.  1,  the  larger  feedback  gains  or  lead  compensation 
necessary  to  improve  the  standard  AOA  system  performance  result  in 
excessively  overdriving  the  engine.   The  longitudinal  feedback  provided 
an  increase  in  performance  and  a  more  efficient  use  of  the  engine.   The 
thrust  response  induced  was  directly  related  to  the  damped  natural 
frequency  of  the  phugoid  mode.   As  the  frequency  was  increased,  as  in 
Figure  12  for  increasing  K   ,  peak  to  peak  thrust  values  increased  as 
well.   However,  for  configurations  in  which  the  damped  natural  frequency 
was  not  appreciably  increased,  as  in  Figures  11  and  14,  peak  to  peak 
thrust  values  were  less  than  those  for  the  standard  configuration. 

B.   EFFECTS  OF  ADDITIONAL  APCS  INPUTS 

As  mentioned  in  Section  III.B,  the  theoretical  analysis  neglected 

AFCS/ACLS  loop  closures.   It  was  shown  that  these  closures  do  not 

appreciably  effect  airspeed  response  for  the  system  under  study.   Phugoid 

damping  coefficients  and  damped  natural  frequencies  were  taken  from 

the  root  locus  plots  of  Figures  11  through  14.   Measurements  of  the 

period  and  the  time  to  one  half  initial  amplitude  of  the  airspeed 

responses  were  then  made  for  the  closed  loop  simulations,  which  included 

AFCS  and  ACLS  closures.   The  damped  frequencies  and  damping  coefficients 

calculated  from  these  responses  were  within  a  few  per  cent  of  those 

predicted  by  the  root  locus  plots.   Thus,  it  was  concluded  that  airspeed 

response  as  controlled  by  the  APCS  was  not  significantly  altered  by 

AFCS  and  ACLS  loop  closures. 

1.   Longitudinal  Acceleration 

Only  a  limited  increase  in  performance  was  achieved  by  the 

implementation  of  u  feedback,  as  only  a  minimal  increase  in  phugoid  damp- 

ing  could  be  achieved. 
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2»   Longitudinal  Velocity 

The  addition  of  u  feedback  increased  the  absolute  phugoid 
damping  but  provided  only  a  minor  increase  in  the  damping  coefficient 
and,  therefore,  no  decrease  in  per  cent  overshoot.   The  most  important 
factor  was  the  increased  damped  natural  frequency.   It  was  considered 
necessary  to  keep  the  damped  frequency  near  the  level  of  that  of  the 
standard  configuration,  as  a  much  higher  value  would  prove  unacceptable 
due  to  the  rapid  oscillations  in  thrust  which  would  result. 

3.  Proportional  Plus  Integral  Longitudinal  Velocity 

The  inclusion  of  an  integral  term  had  an  adverse  effect  on 
airspeed  response  for  an  initial  airspeed  error.   This  was  predicted  by 
the  root  locus  plot  as  shown  in  Figure  13.   The  response  to  a  longitud- 
inal gust  field  was  observed  to  be  improved,  as  documented  in  Table  VII< 
This  shows  up  in  the  root  locus  plots  of  Figures  10  and  13  as  an  approx- 
imate cancellation  of  the  phugoid  poles  with  u/u   zeroes. 

4 .  Longitudinal  Velocity  Plus  Accelc ration 

The  inclusion  of  a  combination  of  proportional  plus  derivative 
of  longitudinal  velocity  in  the  APCS  provided  the  best  increase  in  over- 
all system  performance.   For  values  of  K  greater  than  0.5,  the  damped 
frequency  increases  substantially  with  the  feedback  gain,  K*  .   For  a 
value  of  K  equal  to  0.5  a  large  increase  in  damping  was  achieved  while 
the  damped  frequency  was  held  nearly  constant.   Improvement  of  APCS 
performance  rapidly  diminished  as  K  was  reduced  below  0.5. 

q 

A  stable  first-order  pole  was  shown  to  approach  the  origin  as 
longitudinal  feedback  gains  were  increased  [Figures  12  and  14],  and  its 
presence  was  evident  in  the  digital  simulations.   For  the  configuration 
described  in  Section  IV.B4  with  an  initial  airspeed  error,  a  residual 
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airspeed  error  of  opposite  sign  and  with  a  value  of  ten  per  cent  of  the 
initial  airspeed  error  remained  and  slowly  decayed  after  the  phugoid 
had  essentially  damped  out.   The  effect  of  the  pole  was  considered  to 
be  of  minor  significance.   As  shown  in  Figure  15,  there  was  a  correspond- 
ing zero  in  the  w/6   numerator.   Consequently,  the  first  order  pole  did 
not  appear  in  the  angle  of  attack  response. 

C.   FINAL  CONFIGURATION 

The  choice  of  the  final,  revised  configuration  was  based  on  several 

factors:   phugoid  damping,  damped  natural  frequency  and  gust  response. 

A  combination  of  longitudinal  velocity  and  acceleration  inputs  was  the 

only  combination  able  to  provide  a  substantial  increase  in  phugoid 

damping  while  maintaining  a  constant  damped  natural  frequency.   The 

choice  of  the  value  of  K*  was  somewhat  arbitrary.   For  a  value  of  2.0 

u 

the  effective  longitudinal  velocity  feedback  gain  was  1.0,  which  was  of 

the  same  magnitude  as  K   in  units  of  volts  per  ft/sec;  at  a  higher 

value  of  K*  the  airspeed  input  becomes  dominant  over  AOA.   Improved 

performance  at  a  larger  value  of  K'  was  indicated  by  the  root  locus 

u 

plot  and  was  verified  by  simulation  [Table  VII], 

Longitudinal  gust  response  was  substantially  improved  for  all 
configurations  with  longitudinal  feedback.   The  most  significant  improve- 
ments, as  indicated  by  the  RMS  velocity  errors,  were  achieved  for  the 
configurations  in  which  the  damped  natural  frequency  was  increased  over 
that  of  the  standard  configuration.   The  ratio  of  u  to  u  feedback  chosen 
represented  a  tradeoff  between  gust  response  and  frequency  criteria. 

Vertical  gust  response  as  measured  in  simulation  was  not  appreciably 
altered.   As  that  was  expected  prior  to  the  study  (Section  I.C),  there 
was  no  theoretical  analysis  attempted  in  that  area. 
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For  an  initial  altitude  error  for  all  configurations,  the  induced 
airspeed  error  was  primarily  an  AOA,  or  w(t),  transient  which  decays  in 
five  seconds.   The  induced  airspeed  error,  u(t),  was  much  smaller  in 
comparison.   Thus,  the  addition  of  the  longitudinal  feedbacks  have  a 
negligible  effect  on  responses  induced  by  an  initial  altitude  error. 
Based  on  this  result,  it  was  concluded  that  airspeed  response  at  push- 
over from  level  flight  to  glideslope  acquisition  would  be  unaltered. 

The  frequency  response  plots  of  Figures  19  through  21  indicate  the 

frequency  range  which  was  modified  by  the  APCS  revision.   The  amplitude 

plots  for  u/6   and  u/u   indicate  this  range  was  from  0.07  to  0.5  rad/sec, 
eg 

while  only  the  amplitudes  at  low  frequencies  were  attenuated  in  the 
w/6   plot.   The  additional  favorable  effect  of  an  approximate  three  dB 
attenuation  at  the  phugoid  break  frequency  of  the  revised  configuration 
was  not  indicated  on  the  asymptotic  plots. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.   APCS  CONFIGURATION 

Substantial  improvement  in  performance  is  possible  through  the  use 
of  longitudinal  feedback.   However,  only  an  acceleration  term  can  be 
practically  obtained  with  the  current  APCS,  as  it  would  entail  only  the 
implementation  of  an  accelerometer  signal.   Once  an  airspeed  feedback 
is  considered,  the  need  for  an  airspeed  reference  is  required,  necessitat- 
ing a  major  APCS  redesign.   An  airspeed  reference,  manually  set  by  the 
pilot,  would  not  be  acceptable.   In  addition  to  the  added  pilot  respon- 
sibility mentioned  in  Section  I,  the  combined  AOA/airspeed  system 
references  would  require  airspeed  reference  accuracy  greater  than  the 
pilot  could  be  expected  to  provide. 

In  summary,  only  limited  improvement  in  performance  is  possible 
within  current  limitations.   It  is  recommended  that  the  current  APCS 
design  concept  be  considered  for  revision.   Present  day  aircraft,  such 
as  the  A-7E  and  the  F-14A/B,  have  precise  airspeed  measuring  equipment 
(air  data  computers)  and  digital  computers  on  board.   The  possibility 
of  utilizing  those  devices  to  accurately  measure  and  reference  airspeed 
for  its  incorporation  into  the  APCS  should  be  studied. 

The  possible  improvement  of  gust  response  gained  by  the  revised 
APCS  configuration  is  partially  governed  by  the  method  of  sensing  air- 
speed.  If  an  inertial  device  is  utilized,  the  actual  airspeed,  u,  is 
measured,  but  the  term  (u-u  )  is  measured  if  an  air  data  computer  is 
utilized.   The  effect  of  including  the  gust  velocity  in  the  sensed 
airspeed  signal  was  to  degrade  performance  below  that  gained  by  the 
revised  configuration  (e.g.,  the  RMS  value  of  velocity  error  for  a 
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horizontal  gust  input  was  4.1  ft/sec.     With  the  gust  velocity  included 
in  the  feedback  signal,  the  error  was  5.0  ft/sec,  still  substantially 
improved  over  the  original  value  of  6.3  ft/sec). 

Angle  of  attack  response  was  not  degraded  by  the  inclusion  of 
longitudinal  feedback  but,  in  fact,  was  improved. 

B.  ENGINE  RESPONSE 

The  effect  of  the  nonlinear  engine  time  constant  on  APCS  performance 
was  shown, to  be  minor  for  the  range  of  operating  conditions  considered. 
The  limitations  imposed  on  APCS  performance  by  a  large  engine  lag 
are  much  less  severe  for  a  system  incorporating  longitudinal  velocity 
and  acceleration  feedback  than  for  the  standard  angle  of  attack  config- 
uration. 

C.  THEORETICAL  ANALYSIS 

Restricting  the  multiple  loop  control  system  analysis  technique  to 
the  APCS  loop  closures  was  shown  to  be  valid  for  an  evaluation  of  air- 
speed control.  The  implementation  of  the  method  provided  a  systematic 
approach  for  the  study  of  APCS  performance. 

D.  GENERALIZATION  OF  RESULTS 

The  choice  of  a  single  APCS/aircraf t  system  for  the  analysis  may 
restrict  the  validity  of  some  of  the  results  to  that  particular  system. 
However,  as  all  current  Navy  APCS's  use  the  same  basic  AOA  system,  it 
is  postulated  that  the  general  concept  presented  is  applicable  to  other 
aircraft.   Correlation  of  results  with  other  systems  should  be  attempted. 
The  study  was  restricted  to  the  fully  automated  carrier  landing 
mode.   The  validity  of  the  results  for  a  pilot-controlled  attitude  loop 
should  be  investigated. 
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TABLE   I 
AIRCRAFT  GEOMETRY  AND  STABILITY  DERIVATIVES 


STABILITY 
DERIVATIVES 


X 


u 


X 


w 


w 


w 


M 


M« 


w 


M 


w 


M 


M, 


X 


AT 


AT 


M 


AT 


-  0.05435 
0.064327 

-  0.286953 

0 

-  0.528871 

-  0.000165 

-  0.000289 

-  0.007964 

-  0.327532 
0.732836 

-14.713536 

-  2.188878 

0.001317 

-  0.000250 
0.000004 


APPROACH 
PARAMETERS 

s 

375  ft2 

w 

24,000  lb 

I 

yy 

2 

68,000  slug-ft 

MAC 

10.84  ft 

ra 

746  slugs 

h 

Sea  Level 

P 

0.002378  ^H|£ 
ft 

uo 

218  ft/sec 

ao 

12  deg 

e.g. 

28,6  %  MAC 

1 
X 

6.7  ft  (ahead 

of  e.g.) 
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K, 


K; 


K 


n 


K 


K 


K 


R 


x 


TABLE  II 
CONTROL  SYSTEMS  PARAMETERS 


AFCS 


1.9  rad/rad 
1.0  rad/rad-sec 
3.0  deg/g 


ACLS 


1/5.4 
1.0 
0.376 
1.0 
1/15 
2.5 


K, 


w 


K 


a' 


K 


a. 


a 


K 


APCS 


143.2  volts/rad 
2.3   sec 

2.6  volts/unit  a 
0.55  volts/unit-sec 
0.35  sec 
30.0  volts/g 
1.0  sec 
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TABLE   III 
CSMP  FUNCTIONAL  BLOCKS 


General  Form 


Function 


Y  =    INTGRL(IC,X) 
Y(0)    =    IC 

Integrator 


Y  =   REALPL(IC,P,X) 
Y(0)    =    IC 

st 

1        Order  Lag 


Y  =   CMPXPL(IC    ,IC2,P1,P2,X) 

Y(0)    =    IC 
Y(0)    =    IC^ 

Complex   Pole 


Y  = 


/*Xdt  + 


IC 


Laplace  Transform: 


PY  +  Y  =  X 
Laplace  Transform: 


Ps  +  1 


Y  +  2P  P£Y  +  P2Y  =  X 

Laplace  Transform: 
1 


2  2 

s  +  2P  P2s  +  P2 


Y  =  AFGEN(FUNCT,X) 

Arbitrary  Function  Gener- 
ator(Linear  Interpolation) 


Y  =  FUNCT(X) 


Y  =  RNDGEN(P) 

P  -  ANY  ODD  INTEGER 

NOISE (RANDOM  NUMBER)  GENER- 
ATOR WITH  A  UNIFORM 
DISTRIBUTION 


UNIFORM  DISTRIBUTION  OF 
VARIABLE  Y 

p(Y)  =  Probability  distri' 
bution  Function: 

k 


i.ota 


1.0 


-$*~  Y 
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TABLE   IV 
AIRCRAFT   STATE   EQUATIONS  MATRICES 

[A] 


X  X  -g 

u  w 


Z  Z 

u  w 


0  IL 

1-Z-  1-Z-  0 

w  w 


1.0 


M-Z  M-Z  M«Un 

M+JLH  M  +  wji  o                      M  +.-==? 

u   1-Z*  w   1-Z*  q    1-Z» 

w  w  w 


CB] 


xe  x._  -X  -X 

■6  AT  u  w 

e 


Zf  ZAm  -z  -z 

6  AT  •  u  w 

e 


MwZ6  M-Z  m  M-Z  M-Z 

e  w  AT  w  u  ww 

M     + -  M     +            ■  -M  -  — —  -M  - 

5        1-Z-  AT     1-Z-  u    1-Z-  w   1-Z- 

e            w  w  w  w 
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TABLE  IV  (continued) 


1-Z- 


w 


-1     M  + 


M  Z 
w  u 


x      u      1-Z' 


w 


[c] 

-1.0 


w 


M'Z 
w  w 


-1      M  +  , 
x     w   1-Z' 


w 


w   0 


-1     M  +. 
xl    q    1-Z 


w 


[D] 


M'Z. 

w  o 

e  w 


JAT 


M'  ZA  „, 
MV    1-Z- 


-Z 


M-Z 

1    IM  +f2Jf 

x\    u    1-Z' 
w 


-Z 


w 


M-Z 

V7    W 


1      M  +. 
x\    w    1-Z 


w 


[E] 


0 
1.0 
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TABLE  V 

AIRCRAFT  OPEN  LOOP  TRANSFER  FUNCTION 
NUMERATORS  AND  DENOMINATOR 


s4  +  0.973939s3  +  l.'97797s2  +  0.109846s  +  0.070689 


NAT    =  "  °'00025s  +  0.0003986s   -  0. 0001281s  +  0.00003824 

N6     =   0. 00000407s2  +  0.00000422s  +  0.00000320s 
AT 


u 


0.001317s3  +  0.001195s2  +  0.002435s  -  0.0001321 
AT 


NY     ="  14.7135s3  -  483.007s2  -  26.3803s  -  20.1217 
6e 


N"         0.732836s3  -  0.272704s2  +  40.6515s  +  33.4610 
6 


.u 
e 


NU     =   0.054534s  +  0.0685976s  +  0.109846s  +  0.070689 
u 
g 

N5  AT  =   °-0191945s  +  0.635365s  -  0.0194913 
e 

N"  ®   =   0.00288013s  +  0.0013343 
e 

w    u  2 

Ne    l\n     =  -  0.0191945s      -  0.635365s  +  0.0194913 
6   AT 
e 

nY  ?m     =   -  0.00060607s    -   0.00085794 
6  AT 
e 

NU  "        =   -   0.00039155s      -   0.00012806s  +  0.00003824 
u  AT 
g 

NU   ,9m     =        0.0000001131s  +  0.00000320 
u  AT 
g 
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TABLE  VI 

AIRCRAFT  TRANSFER   FUNCTION  NUMERATORS  AND 

DENOMINATOR   FOR   STANDARD  APCS   LOOP  CLOSURES 

T     =1.8    SECONDS 
e 

A"      =      S(s+l)(s+2l857)(s+0.56){^0-303'0'278^-0-263) 
(0.032) (-1.099) [0.346, 1.44] (-2.89)1 


6  s(s+l)(s+2.857)(s+0.435)(s+0.56)  \ £°-504>°- 1823 

(0.150) [-0.176, 1.38] (-0.807) (-3. 41) (49. 5)1 


N 
std 


<]    =     rsTT)T^o!;35)(S+0.56){^-389'0-257^-0-253> 

std  v 

(-0.584) (-1.05) (-32.2)V 


Ul    =        ,  ,^s0/l^5l^s,    ,n'V^s  <(  [0 .  308 , 0 . 7 16]  (0 .  00948) 
u  s(s+l) (s+2.857) (s+0.56)    (  JV  ' 

[0.0836, 1.74] (-1.58) (-3.37)1 


N 

1 

std 


2        2 
Notation:  [f    ,  co    represents  (s  +2fu)  -hu  ) 
^    n  s  n  n 


(t)       represents  (s-t) 

positive  £  and  negative  T  represent  stable  roots 
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TABLE  VII 
RMS  VELOCITY  ERROR  SUMMARY 


APCS 
Configuration 

"o=  I 

wQ=-5 

wQ-  5 

Z  =10 

eo 

Z  =-10 

eo 

w  =5 
g 

u  =5 
g 

Standard 

3.09 

2.87 

3.14 

3.06 

7.40 

6.26 

Long  Accel 
K'=1.0 
"=2.0 
=3.0 

3.02 
3.03 
3.09 

2.85 
2.91 
3.00 

3.08 
3.05 
3.03 

3.03 
3.01 
2.96 

7.01 
6.75 
6.57 

5.63 
5.28 
5.05 

Long  Veloc 
K  =1.0' 
U=2.0 
=3.0 

2.78 
2.69 
2.59 

2.57 
2.38 
2.20 

3.09 
3.05 
3.01 

3.05 
3.03 
3.01 

7.11 
6.72 
6.40 

4.16 
3.44 
3.28 

Prop+Int  u 
K  (u+iffu) 
nz.   =1.0 

"=3.0 

3.26 

3.88 

2.95 
3.05 

3.16 
3.10 

3.11 

3.11 

7.79 
6.23 

4.46 
3.30 

Prop+Accel  u 

K*(u+K  u) 
u    q 

K  =1.0 

V=1.0 

"=3.0 

K  =0.5 
V-1.0 

"=2.0* 
=3.0 

K  =0.2 
qK-=1.0 
"=3.0 

2.65 
2.35 

2.75 
2.61 
2.52 

2.88 
2.73 

2.48 
2;30 

2.60 
2.48 
2.41 

2.72 
2.64 

3.05 
2.99 

3.07 
3.03 
3.01 

3.07 
3.02 

3,02 
2.98 

3.02 
3.00 
2.96 

3.03 
2.99 

6.76 
6.18 

6.87 
6.55 
6.35 

6.94 
6.46 

4.07 
3.55 

4.61 
4.06 
3.84 

5.12 
4.33 

*  Underlined  values  denote  revised  configuration. 
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FIGURE  2.    TF41-A-2  ENGINE  CHARACTERISTICS 
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AFCS   EQUATION 
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FIGURE  4.   AFCS  AND  ACLS  TRANSFER  FUNCTIONS 
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APPENDIX  B 
GUST  MODEL  ANALYSIS 

To  simulate  a  gust  input  for  the  CSMP  simulation,  a  random  number 
generator  was  used.   The  CSMP  function  RNDGEN  has  a  uniform  probability 
distribution  function  with  possible  values  from  0.0  to  1.0.   As  the 
desired  mean  of  the  gust  input  was  0.0,  a  constant  value  of  0.5  was 
subtracted  from  the  RNDGEN  output.   The  output  of  the  random  number 
generator  was  passed  through  a  CSMP  second-order  filter  with  a  break 
frequency  of  1.0  rad/sec  and  a  damping  coefficient  of  0.707  and  was 
multiplied  by  a  constant  to  achieve  an  RMS  value  of  5.0.   The  following 
is  a  power  spectral  analysis  of  the  gust  model. 

Denoting  the  RNDGEN  output  signal  in  the  frequency  domain  as  X(ju)), 
the  filter  Fourier  Transform  as  H( jau)  ,  and  the  resultant  gust  signal  as 
Y( joo)  ,  the  relationship  between  X(jo>)  and  Y(jcu)  is  shown  in  Figure  B.l. 
The  output  power  spectral 


X(ju>) 


^  Y(juj) 


FIGURE' B.l   GUST  MODEL  BLOCK  DIAGRAM 

density  is  equal  to  the  input  power  spectral  density  multiplied  by  the 
absolute  value  squared  of  the  filter  Fourier  Transform: 


$  ( jco) 

yy 


H(ju>) 


XX  J 


(1) 


The  mean  square  value  of  the  output,  or  average  output  power,  can  be 
determined  by  evaluating  the  autocorrelation  function,  cp   (t)  ,  at  j   =   0: 

yy 
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y  "  -  cp   (0)  =  jjj-  J°V  (jw)  dco 


=  5F  r|H(jU))f  ^(J")*" 


(2) 


y  is  defined  as  the  RMS  value  of  y(t). 


A  typical  time  plot  of  x(t)  is  shown  in  Figure  B.2. 


x(t) 


>  t 


FIGURE  B.  2   RNDGEN  Output  vs.  Time 

The  random  number  generator  output  is  a  discrete  value  at  each  At.   The 
resultant  curve  through  the  points  represents  the  fourth-order  Runge- 
Kutta  approximation  of  the  function.   By  approximating  the  curve  by  a 
series  of  stops,  a  rough  approximation  of  the  autocorrelation  function, 
cp   (t)  ,  was  determined. 

The  autocorrelation  function  is  defined  by  the  equation, 


cp   (T)   =  f   x(t)  x(t  +  t)  dl 


(3) 


For  the  step  series  function  utilized  to  approximate  the  RNDGEN  output, 
each  segment,  At  in  length,  is  uncorrelated  with  every  other  segment; 
that  is. 


CPXX(T)  =  0 


|  Tj  >  At 


(4) 


To  determine  cp      (T)    subject   to   the   restriction  of     Equation    (4),    x(t) 


xx 


was  depicted   as   in   Figure   B.3   [Ref.    B.l].      By   applying   Equation    (3) 


A   x(t) 
a 


-i >   t 

At 


At 
2 


FIGURE  B.3   Plot  of  x(t)  to  Determine  cp   (t) 

XX 


to  Figure  B.3  and  requiring  that  cp   (t)  evaluated  7  =  0  be  equal  to  the 

XX 

variance,  cp   (t)  was  determined. 

XX 


CO       (T)     =   o 

Yxx  x 


2/_  Jrl 

At 


+  1 


<     At 


(5) 


The    power   spectral   density   can  be   expressed 


*      (j«0      =      f%      (T)e"jUjT   dr 
XX  J        Txx 


(6) 


By  substituting  Equation  (5)  into  Equation  (6)  and  noting  that  cp   (T) 

XX 

is   even   function, 


2  1 

<^x(jw)    =  2ax    (-     —  T  +   1)    coscoT 


(7) 


2a 


2, 

co  At 


(    1    -   coscuT) 


By   successively   applying  L'Hospital's   Rule 


E>     (0)      =     a     At 

XX  X 


(8) 
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<3?  (jco)  is  periodic  with  a  period  equal  to  2rr/At.   For  a  small  integra- 
tion  step  size  with  the  filter  cutoff  frequency  of  1.0  rad/sec,  $  (jco) 

XX 

is  essentially  flat  in  the  frequency  range  of  interest  and  can  be 
considered  a  constant  at  the  zero-frequency  value  shown  in  Equation  (8) . 
For  a  second-order  filter  with  a  break  frequency  of  1.0  rad/sec,  the 
magnitude  of  the  transfer  function,  H( jcu) ,  squared  has  a  -40  dB/decade 
slope  for  frequencies  greater  than  1.0  rad/sec.   Based  on  the  flat  power 
spectral  density  at  low  frequencies  and  the  filter  cutoff  frequency, 
the  evaluation  of  Equation  (2)  was  approximated  by  taking  <£  (jco)  as  a 

XX 

constant  at  its  zero-frequency  value.   The  resulting  integral  was 

2 

2    1    °°  m 

V<0)   =  CT*  ''^1      ~  ,2    —  -  ; —  *°  w 

-00   (jco)   +  2lo)    (jco)  +  co 
n  J      n 


and  was  of  a  form  tabulated  in  Ref.  B.2.   Substituting  for  to  and  K   and 

n 

evaluating, 

cp      (0)      =     0.354ct  2«   At   =  y  2 
Yyy     '  x  J 


At  y  =  5.0   ft/sec   and  At   =  0.05   sec, 


a     =  37.6   ft/sec  (10) 

x 


In  actual  simulation  the  value  of  a     was  strongly  dependent  on  the 

RNDGEN  used  in  the  simulation,  varying  from  1.7  to  3.1  with  a  mean  of 

2.4,  all  values  falling  below  the  predicted.   RMS  values  of  the  input  to 

the  filter  (the  outputs  of  the  RNDGEN' s)  were  essentially  constant  at 

the  predicted  value.   The  underestimation  of  the  value  of  a     required 

x 

to  achieve  the  desired  value  of  a  was  a  direct  result  of  the  step 

y 

approximation  to  the  fourth-order  integration.   The  reason  for  the  large 
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variance  in  output  RMS  values  was  not  determined  but  was  suspected  to 
be  caused  by  the  short  run  time.   The  length  of  each  run  was  30  seconds, 
or  3600  integration  points.   In  the  final  model,  a     was  adjusted  to 
obtain  the  desired  value  of  5.0  ft/sec   for  <j   . 

y 
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APPENDIX  C 
MULTIPLE  LOOP  CONTROL  SYSTEM  ANALYSIS 


The  multiple  loop  analysis  technique  described  below  is  derived 
and  demonstrated  for  a  general  case  in  Ref.  C.l.   The  technique  is  a 
procedure  for  expressing  specific  closed  loop  transfer  functions  of  a 
multiple  loop  system  in  terms  of  the  elements  of  the  open  loop,  feed- 
forward and  feedback  matrices.   The  following  is  a  two  control  input, 
three  output  example. 

Consider  the  system 


[A] 


{X} 


[B] 


{6}   +     [C]       {R} 


Hi'    312'    ai3 

(  Xx(s) 

a21'    a22'    a23 

X2(s) 

a31'    a32'    a33 

(x3(s) 

bll>  bl2 


h2V   b22 


b31,  b32 


6^3) 


i2(s)j  + 


Cll'  °12 


C21'  °22 


c31'  °32 


R^s) 


R2(s) 


(1) 


Matrix  elements  are  general  and  may  contain  operators.   [A]  is  the  open 

loop  system  matrix,  {x}  is  a  vector  of  variables,  {6}  is  a  vector  of 

control  variables,  and  [r]  is  a  vector  of  external  disturbances.  The 

system  will  contain  one  feedforward  and  two  feedback  paths  to  drive  the 

control  variable  5  ',  as  depicted  in  Figure  C.l. 

Rl    *2 
J^ 

>  xx 

>  xn 


>  X, 


FIGURE  C.l      Control    System   Block  Diagram 
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The  control  equation  is 


61  =  -hlX2  -  h2x3  +  Yl62 


(2) 


The  (s)  notation  has  been  dropped  for  simplicity.   Substituting  Equation 
(2)  into  (1)  and  rearranging 


an,  a12+h1bn,  a13+h2bu 


a21,  ^2+l^b^,  a23+h2b21 


a31,  a32+h1b31,  a33+h2b31 


X, 


b12+Y1bu 


b22+Ylb21 


b32+Ylb31 


52  + 


Cll'  C12 


C21'  °22 


c     c 

c31'   32 


R, 


(3) 


Equation  (3)  is  represented  by  the  matrix  notation 

[A"]  {X}   =   [B1]  {62}   +  [C]  [R] 

where  A1 '  denotes  two  loop  closures  and  B1  denotes  one  feedforward  path. 
The  closed  loop  system  determinant  can  be  expanded  into  a  sum  of  four 
determinants 


=  A"  = 


all'  al2'  al3 


a21'  a22'  a23 


a31'  a32'  a33 


all'  a12'  h2bll 
a2l,  a22,  h2b21 


a31'  a32'  h2b3l 


all'  hlbll'  a13 
a21'  hlb21'  a23 
a31'  hlb31'  a33 


+ 


(4) 


au,  h1bn,  h2bn 

a21>  hlb21>  h2b21 
a31,  h^^  h2b31 
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The  determinant  with  a  subscript  of  1  is  the  open  loop  system  determinant, 
Number  2  is  the  numerator  of  the  open  loop  X9/6,  transfer  function 
multiplied  by  h  .   Number  3  is  the  numerator  of  the  open  loop  X  /6 
transfer  function  multiplied  by  h„.   Number  4  is  identically  equal  to 
zero,  as  the  second  and  third  columns  are  proportional. 

Nc  denotes  the  numerator  of  the  open  loop  X/6(s)  transfer  function. 
Using  this  notation }   Equation  (4)  becomes 


A"  =  A  +  h^  +  ty^ 


(5) 


Equation  (5)  can  be  expanded  for  the  general  case  to  any  order  and 
number  of  closures  as 


An+m     A 


n  -— „  m 


V   V       x 

2-   2.    h..N?j 


(6) 


h..  represents  the  feedback  transfer  function  from  6.  to  x.. 


< 


x„ >  x „  *  6^ 


is  defined  as  the  closed  loop  numerator  of  the 


X./69  transfer  function  for  the  above  loop  closures. 
From  Equation  (3) 


■? 


x2,x3,82  ~*  Sl 


b12  +  Ylbll'  al2  +  hlbll'  a13  +  h2bll 


b22  +  Ylb2l'  a22  +  hlb21'  a23  +  h2b21 


b32  +  Ylb31'  a32  +  hlb31'  a33  +  h2b31 


(7) 


85 


Expanding  Equation  (7) 


N. 


x2 '  x3 '   2     1 


bl2'  al2'  a13 
b22,  a22,  a23 

b32,  a32,  a33 


+ 


Ylbll'  a12'  a13 
Ylb21'  a22'  a23 
Ylb31'  a32>  a33 


+ 


bl2,  h1bn,  a13 

b22'  hlb21'  a23 
b32'  hlb31'  a33 


+ 


b12,  a12,  h2bn 

b22'  a22'  h2b2l 
b32,  a32,  h2b31 


+  0  +  0  +  0 


(8) 


Inspection  of  Equation  (7)  reveals  that  three  determinants  have  dropped  out. 
Equation  (8)  can  be  rewritten  as 


X ..  X  _  X_X/>  X  -  X  ^ 

=  Ne   +  Y.N.   +  h.N  \   +  h.N  e 
62         °1  2  1  2  1 


X2 '  X3 '   ? 


(9) 


where 


xxx2 
6262 


bll'  bl2'  al3 
b21'  b22'  a23 
b3T  b32'  a33 


Again,  Equation  (9)  applies  in  general. 

In  the  following  case  the  output  variable  is  X9,  one  of  the  feedback 
variables. 


N. 


;] 


xo  »  X3 '   9 


aiV    b12+Ylbll>  ai3+h2bll 
a21,  b22+Y1b21,  a23+h2b21 

a31'  b32+Ylb3l'  a33+h2b31 
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In  the  above  notation 


N. 


x2'  x3'  62  ~*  61 


(10) 


In  Equation  (10)  another  crossproduct  was  lost.   This  was  a  special  case, 

but  the  notation  of  Equation  (9)  applies,  noting  that  the  lost  term 

x~x_ 
would  have  been  h.N.  .  .   A  numerator  of  this  form  makes  no  sense  and  is 

defined  as  identically  zero. 

The  following  is  a  numerator  of  a  transfer  function  of  an  output 

variable  due  to  an  external  disturbance. 


::] 


X9  '  X"3  '   9 


c.u,  al2+h1bn,  a13+h2bu 
c2l,  a22+h1b21,  a23+h2b21 
c31,  a^+h^,  a33+h2b31 


X-  X-X«  1   x 

r.      1  r  6,      2  r  6 


(11) 


Equation  (11)  is  similar  to  Equation  (9),  but  the  feedforward  term  is 

missing.   The  notation  of  equation  (11)  applies  to  any  external  input 

for  which  there  is  no  feedforward  term,  either  for  another  control 

variable  or  for  an  actual  external  disturbance,  as  illustrated. 

It  should  be  noted  that  N     and  N     denote  the  same  determinant; 

'  6162      5261 
i.e.,  rows  are  not  interchanged.   The  latter  notation  would  arise  if 

X3l 

N.  were  formed. 

62j 

X2  '  X3  '   2  "" *   1 
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